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$_l , We study the effects of the interaction between the Chiral condensate and the Polyakov loop 

0- on the chiral transition within an effective Lagrangian. We find that the effects of the interaction 

' change the order of the phase transition when the explicit breaking of the Zn symmetry of the 

Polyakov loop is large. Our results suggest that the chiral transition in 2-flavor QCD may be first 

order. 
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I. INTRODUCTION 



o , 

QCD matter with Nf flavors, N colors and zero baryon chemical potential undergoes two finite 
temperature phase transitions, the chiral transition and the deconfinement. The chiral transition is 
associated with the spontaneous breaking of the chiral symmetry, SU(Nf) l x SU (Nf)n — > SU(Nf)y, 
below the critical temperature (T x ) for massless quarks. The order parameter for this transition 
' is the chiral condensate. On the other hand the deconfinement transition is associated with the 

Oh, spontaneous symmetry breaking, Zjy C SU(N) —> 1, above the critical temperature Td for infinitely 

heavy quarks. The order parameter for this transition is the Polyakov loop expectation value. For 
finite non-zero quark masses both the chiral symmetry and the Zn symmetry are explicitly broken. 
• Nevertheless these two transitions show up as crossover, second or first order transitions depending 

on the values of the quark masses. The chiral transition depends on the number of quark flavors 
and the deconfinement transition depends on the number of colors. Since the chiral symmetry and 
the Zn symmetry seem nothing to do with each other one would expect that these two transitions 
occur independently. However lattice QCD calculations have shown that both these transitions occur 
simultaneously, T x — Td [1-5]. Furthermore strong correlation between the chiral condensate and 
the Polyakov loop is observed around the phase transition point This is clear evidence that there 
is interaction between these two order parameters. So studying the interaction between the chiral 
condensate and the Polyakov loop is fundamental to understanding the interplay between the chiral 
transition and deconfinement. There are several studies on the possible causes of the simultaneous 
occurrence of the chiral transition and the deconfinement transition. Mixing between the gauge 
and matter field operators has been suggested to explain the simultaneous chiral and deconfinement 
transitions @>I3- Some other studies consider that for small quark masses the chiral transition drives 
the deconfinement transition 0,0- 

Though most of the studies are concerned with why T x = Td only a few have considered the effect 
of the interaction on the phase transition itself [Tcl | . It seems natural to expect that if the interaction 
between the two order parameters can result in the simultaneous occurrence of the two transitions 
then the interaction may also have important effect on the phase transition. One of the most 
interesting cases to study for the possible effects of this interaction is the 2-flavor chiral transition. 
Lattice QCD calculations have not yet been able to settle on the order of this phase transition. 
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Lattice calculations by different groups do not agree on the order of this phase transition. Some 
lattice groups find the transition is second order Q and other groups find the transition is first 
order [ijj ■ Conventionally this transition is believed to be second order and in the universality class 
of 0(4) Heisenberg magnet 12] . But the effect of interaction between the chiral order parameter 
(<&) and the Polyakov loop (L) may change the behavior of this transition. So in the present work 
we investigate the effect of interaction between $ and L on the 2-flavor chiral transition within an 
effective Lagrangian. 

Previously, the effect of interaction between the chiral order parameter and the Polyakov loop 
has been studied in the renormalization group approach 0. The main difference between the 
present work and previous studies is that we consider the explicit breaking of the Zn symmetry. 
Explicit breaking of the Zn symmetry can be introduced in the effective Lagrangian by terms such 
as ~ (L + L'), ~ (L + L')$f$. Previous studies [l(| have considered interaction term, such as 
(LL^)(&&), which respects both the chiral and the Zn symmetry. However, in the chiral limit, the 
interaction terms need not respect the Zn symmetry. So terms such as ~ (L + L^)$t<j) should be 
considered. As we discuss later such a interaction term is always present if the explicit Zn symmetry 
breaking is large, for example in the chiral limit. For simplicity we consider N = 2 color QCD. We 
expect that different TV will not qualitatively change the physics we are discussing here. For Nf = 2 
the chiral order parameter $ is a four component vector field whereas for N = 2 the Polyakov loop 
L is a real scalar field. In this work we basically study the effect of the three terms, 
and L 2 <I>^<E> in the effective Lagrangian. Our main result is that the strong explicit breaking of the 
Zn symmetry can make the chiral transition first order. We show that the chiral transition can 
be first order even at the mean field level. We also carry out numerical Monte Carlo simulations 
which confirm the first order phase transition for large enough explicit Zn symmetry breaking. We 
mention here that the Nf = 2 chiral transition can be first order from the interaction term L 2 $^$ 
without Zn symmetry breaking [13J. However lattice QCD results indicate that effect of the Zn 
symmetric interaction term is small. A possible first order chiral transition can result more likely 
from the explicit Zn symmetry breaking as we will argue later. 

Conventionally explicit symmetry breaking weakens a phase transition. But our results suggest 
that for a system with two order parameter fields explicit symmetry breaking can make the transition 
stronger. It is interesting to note that the chiral order parameter $ does not couple to gauge fields 
directly. The gauge fields seem to affect the chiral phase transition indirectly, through the Polyakov 
loop. We mention here that the effect of explicit symmetry breaking discussed here should not be 
restricted to the Zn symmetry. We expect that the explicit breaking of chiral symmetry may have 
some effect on the deconfinement transition in the large quark mass region. We mention here that 
interaction between the chiral order parameter and the diquark fields is considered to study the 
chiral/color-superconducting transition at low temperature and high density |14| . 

This paper is organized as follows. In section-II we describe the effective Lagrangian and discuss 
the effect of the interaction between $ and L on the chiral transtition. We describe our numerical 
Monte Carlo calculations and results in section-Ill. The discussions and conclusions are presented 
in section-IV. 



II. THE EFFECTIVE LAGRANGIAN AND THE EFFECT OF THE BROKEN Z N 

SYMMETRY 

We consider the following Lagrangian 9] in 3 dimension for the $ and the L fields, 
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C = i|V$| 2 + i(VL) 2 + F($,L), 

V(*,L) = ^\^ + ^+^ + ^- gL ^- cL ^- eL (1) 

Some of the parameters of this reduced 3D Lagrangian depend explicitly on temperature. This 
Lagrangian is invariant under 0(4) rotation of the <& field. The last two terms of V($, L) break the 
Z2 symmetry (L — ► —L) explicitly. The interactions between the chiral order parameter and the 
Polyakov loop are taken care by 5th and 6th terms in V($>, L). 

The signs of the couplings g and c decides the correlation between the fluctuations of the <E> and 
the L fields. For example when c > 0, the thermal average of the correlation between the $ and L 
fluctuations, ((8L)(6\$\)), is positive. When c < 0, ((5L)(6\$\)) is negative. Such "anti-correlation" 
is seen between the fluctuations of the the chiral condensate and the Polyakov loop in the results 
of lattice QCD calculations 0. The correlation between the variations of the two order parameters 
with respect to temperature is of the same sign as the correlation between the fluctuations. Note 
that in the above Lagrangian when c ^ the <& field acts like an ordering field for the L field. In 
the chiral symmetric phase the chiral order parameter is small and the Polyakov loop expectation 
value is large. The large expectation value of L, in the chiral symmetric phase, can result only from 
the last term in V(<&, L) (with e > 0) because |$| is small. 

As we have mentioned before previous studies |l0j have considered only the interaction term 
g_L 2 |$| 2 . When the coupling parameters c = = e, the chiral transition and the deconfinement 
transition do not always occur simultaneously. For large value of the coupling g the transitions can 
occur simultaneously and are of first order when coefficients of |$| 2 and L 2 are negative in Eq.l 
|l3l IT^ I . A large positive g would increase the critical temperature for the deconfinement transition 
as the coefficcnt of I? term is negative for high temperatures. Lattice results on the other hand 
show that inclusion of dynamical quarks decrease the deconfinement transition temperature. So in 
QCD the coupling g should be small. For small g both the transitions are second order and do not 
occur simultaneously. 

The coupling parameters c and e represent the strength of the Zn explicit breaking. So they 
should increase with decreasing quark masses as the Zn breaking becomes severe. This can be seen 
explicitly in the large quark mass region [H^. To see the effect of the explicit Zn breaking let us 
consider g = = e and c 7^ 0. At the mean field level one can consider the temperature variation of 
the parameters m\ and m 2 L . For simplicity we fix A^,i > 0, m| > and vary the m| parameter. 
To find the m% dependence of L and $ expectation values one has to solve the following coupled 
equations, 

A< i) |$| 3 + m||$|-2c|$|i = 
\ L L 3 + m 2 L L - c\<I>\ 2 = 0. 

We have checked numerically that for large enough c these equations give two solutions which 
correspond to a degenerate minima of the effective potential V($, L) at some particular value of 
m%. It may seem surprising that the potential V($,L) has degenerate minima even though there 
is no cubic term for $ and L in it. However because of the coupling c these two fields are mixed. 
Though the mixing anle varies as m\ is varied. With variation of m\ the minimum of V(<&,L), 
in the |$| — L plane, moves in directions other than the |<&| or L axes. To understand how the 
minimum of V($, L) behaves one should express V($>, L) in terms of variables which are the linear 
combinations of |$| and L. This would invariably lead to cubic terms of the new fields in the effective 
potential. For some choice of parameters the cubic term may then be important to cause degenerate 
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minima of V($, L). Even if the explicit symmetry breaking is not strong enough at mean field level 
fluctuations at higher order can make the transition first order. At one loop the fluctuations of the 
$ field will contribute to a non zero 3-point function of the L field. This 3-point function can be 
calculated perturbatively. In the high temperature approximation the 3-point function is given by, 

c 3 

~T— , (2) 

m A 

assuming zero momentum to all the external L— lines, m here is the mass of the $ field fluctuations. 
Given a suitable choice of parameters in the effective Lagrangian the three point function can be 
significant. The consequence of this is a a(T)L 3 term in the potential V($, L) with temperature 
dependent a(T). This term can cause discontinous change in L as temperature varies. When L 
field changes discontinuously the $ field will also changes discontinuously because of the coupling 
term ci$ 2 . If the coefficient of the $ 2 term changes from +ve to —ve due to discontinuous change 
in L then $ will change discontinuously from zero to non-zero. The L 3 term can also come from 
other types of explicit symmetry breaking coupling terms but we think c_L$ 2 is the simplest term 
in our model. 



The situation with c = and g ^ ^ e is same as the one discussed above. When the explicit 
symmetry breaking parameter e is large, Z 2 symmetry of the L field is lost and the L field always 
has non-zero expectation value L a . In order to study the fluctuations one must expand the potential 
V($,L) around L 0} L = L a + L. This gives rise to a term like gL L& 2 coming from the expansion 
of gL 2 <& 2 around L — L - Now the term gL L$ 2 is similar to the one discussed above. So there 
can be first order transition when g ^ ^ e like in the case when c ^ 0. We observed that at the 
mean field the transition becomes second order for large quartic couplings A$ . Now in the following 
section we discuss the numerical Monte Carlo simulations of the effective Lagrangian (Eq.l) and the 
results. These calculations include higher order as well as non perturbative fluctuations. 



III. NUMERICAL CALCULATIONS AND RESULTS 



The numerical Monte Carlo calculations of the the model Eq.l is done by discretizing the action 
S = J Cd 3 x on a 3 dimensional N 3 lattice. We employ the following rcscaling of the fields variables 
and the couplings, 



<&J~n$ A$ 2 2 — 4A$ — 

<P — ► ,A$ — ► -o-,TO$ — ► 5 ; 

L^/kE A L 2 2 - 4\ L - 6k l 

9 c 



e--^=- (3) 



a is the lattice spacing and k$ and Kl are the hopping parameters for the $ and the L fields 
respectively. After the rescaling of the fields and the coupling parameters the discretized lattice 
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action becomes, 

X 

- kl L x L x+l i + L 2 X + A L (L 2 - l) 2 

9^|$ x | 2 - ci^l^l 2 - ci x (4) 

Here ^(L^) represents the value of the field at the lattice site x. x + [i represents the six 

nearest neighbor lattice sites to x. We adopted the pseudo heat-bath method used for the Higgs 
updating in SU(2)+Higgs studies 0]. To update <& x at a lattice site x we write the probability 
distribution P($ x ) of Q x as, 

P(4 a ) ~ Exp[-Si($ x )-S 2 ($ x )}, with 



A x 2 
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x 1 a cL-r + qL 2 , . 

A = K * ^> + - B = 1 -2A; + 2^ + ^A^ (5) 

The coefficient a is a parameter chosen so that we get a reasonable acceptance rate for the new 
$ x . Once a is chosen, a Gaussian random number is generated according to the distribution 
Exp[— 5i(<& x )]. Then this random number is accepted as the new value of <fr x with the proba- 
bility Exp[— 52 ($ x )]- Using this procedure & x is updated at all the lattice sites. Then we do the 
updating of L x along the same steps. The process of updating <& x and L x on the entire lattice is 
repeated about 20 times between successive measurements. We measure the magnetizations 

M * = ^E $ *> M ^ = ^E L - ( 6 ) 

X X 

where V = Nf is the number of lattice sites. The expectation values of the $ and the L fields are 
given by the thermal averages (average over the measurements), ($) = (|M$|) and (L) — (Ml) . We 
take the absolute value of M$ for (<&) because a normal average of M$ is usually not a well behaved 
observable. 

The numerical simulations were carried out on a N s — 16 lattice. In this work we do not intent to 
explore the phase diagram of the model (Eq.l) but to show that for suitable choice of parameters 
the phase transition can change from second order to first order. Here we present results for two 
sets of parameters. For one set we fix the couplings g = e = and for the second set we fix the 
coupling c = 0. For the first set of parameters we choose, A$ = 0.004, \l — 0.0020, Kl = 0.01 and 
c = 0.1. We observe the hysteresis of ($) and (L) by varying the parameter In FIG.l we show 
the hysteresis loop of ($) and in FIG. 2 we show the hysteresis loop of (L). 

Since we take c to be positive we see (<&) and (L) increase or decrease simultaneously. The choice 
of the values for the parameters is such that the variation of ($) and (L) are similar in magnitude. 
For c < increase in ($) should lead to decrease in (L). So in this case the hysteresis loop for (<f>) 
will look somewhat similar to Fig.l while the hysteresis loop for (L) will be more or less inverted 
about the y-axis. 

For the second set of parameters we choose A$ = 0.0055, Al = 0.0010, kl = 0.14, g = —0.02 and 
e = 0.9. The hysteresis loop of ($) and (L) are observed by again varying the parameter The 
choice of kl and e is such that the expectation value of L is non-zero and positive. The sign of g 
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FIG. 2: The hysteresis of (L) vs k$ 



assures that increase in (<&) leads to decrease in (L) and vice-versa as the parameter k$ is varied. 
The hysteresis curves of the two order parameters are shown in FIG. 3 and FIG. 4. The values of 
A$,l in our calculations are chosen so that we can see first order phase transition clearly and the 
variation of ($) and (L) are of order 0(1) across the transition point. Note that with suitable choice 
of kl and Xl one change the average of the Polyakov loop across the transition point. The choice 
g was such that the chiral transition turned out to be second order when the coupling e was set to 
zero. We also did simulations on a 4D lattice. The results in this case are very similar to the 3D 
simulations. 
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FIG. 3: The hysteresis of ($) vs k$ 



FIG. 4: The hysteresis of (L) vs k$ 



The results in FIG. 1-4 show strong first order phase transition for the $ and L fields. By fixing 
the coupling parameters g, c and e we observed that the strength of the transition depends on the 
values of the quartic couplings A$ and A^. The transition becomes weaker with increase in any 
of the quartic couplings A$.l. However even for larger quartic couplings a suitable choice of the 
parameters <?, c and e made the transition strong first order. 
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We also did calculations with small explicit symmetry breaking for the <& field by considering a 
linear $ term in the Lagrangian. In this case we found that the hysteresis loops of both ($) and (L) 
shrinking in size with increase in the coefficient of linear <E> term in the Lagrangian. These results 
suggest that the transition becomes weaker when the chiral symmetry is explicitly broken. 

IV. DISCUSSION AND CONCLUSIONS 

Using a simple effective Lagrangian, which captures the chiral symmetry and Zn symmetry of 
QCD, we have investigated the effect of the explicit symmetry breaking on the chiral phase transition. 
Since we consider the chiral limit the Zi symmetry of the Polyakov loop is explicitly broken. As we 
incorporate the explicit Z2 breaking interaction terms in the effective Lagrangian we find the chiral 
transition becomes first order. We observed that the first order transition becomes weaker when 
a small explicit symmetry breaking is considered for the chiral order parameter <i>. The results of 
our calculations show that two flavor chiral transition is of first order for large enough Z^ explicit 
breaking. As we have mentioned before at present some lattice studies suggest that the transition 
is first order ^lj an d some other studies show the transition is second order ||. These conflicting 
results may be because the quark masses studied are not small enough or the lattices used are not 
big enough. 

In the previous study for Nf — 2 and N — 3 ^(j > the second order chiral transition becomes first 
order when coupled to the Polyakov loop. This is because the dcconfincmcnt transition is first order 
for N — 3 with a cubic term ~ (L 3 + L* 3 ) in the effective potential with exact Z3 symmetry. However 
when the quark masses are finite and decrease the deconfinement transition becomes weaker. This 
can be understood due to explicit breaking of Z3 symmetry. Already in the large quark mass region 
the deconfinement transition becomes crossover which implies the the explicit symmetry breaking 
dominates over the effects of the above Z 3 symmetric cubic term. For smaller quark masses the 
explicit Z3 symmetry breaking likely grows and expected to be maximal in the chiral limit. So 
for smaller quark masses, i.e in the chiral limit one should rather consider the effect of the explicit 
breaking of the symmetry The effects of explicit symmetry breaking discussed in this work 
should not be restricted to the explicit breaking of the Zm symmetry. For 2-flavor and 2-color 
QCD the deconfinement transition in the large quark mass region may have the effects coming from 
the explicit breaking of the chiral symmetry. It may be possible that this effect change the phase 
transition behavior of the deconfinement transition in the heavy quark mass region. 
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